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ABSTRACT 

We obtained stellar ages and metallicities via spectrum fitting for a sample of 575 bulges with spectra 
available fr om the S l oan Dig ital Sky Survey. The structural properties of the galaxies have been studied 
in detail in lGadottil ()2009bl ) and the sample contains 251 bulges in galaxies with bars. Using the whole 
sample, where galaxy stellar mass distributions for barred and unbarred galaxies are similar, we find 
that bulges in barred and unbarred galaxies occupy similar loci in the age vs. metallicity plane. 
However, the distribution of bulge ages in barred galaxies shows an excess of populations younger 
than 4 Gyr, when compared to bulges in unbarred galaxies. Kolmogorov-Smirnov statistics confirm 
that the age distributions are different with a significance of 99.94%. If we select sub-samples for 
which the bulge stellar mass distributions are similar for barred and unbarred galaxies, this excess 
vanishes for galaxies with bulge mass logAf < 10.1 while for more massive galaxies we find a 
bimodal bulge age distribution for barred galaxies only, corresponding to two normal distributions 
with mean ages of 10.4 and 4.7 Gyr. We also find twice as much AGN among barred galaxies, as 
compared to unbarred galaxies, for low-mass bulges. By combining a large sample of high quality 
data with sophisticated image and spectral analysis, we are able to find evidence that the presence of 
bars affect the mean stellar ages of bulges. This lends strong support to models in which bars trigger 
star formation activity in the centers of galaxies. 

Subject headings: galaxies: active — galaxies: bulges — galaxies: evolution — galaxies: stellar 
content — galaxies: structure 



1. INTRODUCTION 

A large number of studies have investigated the im- 
pa ct of bars on the evolution of galaxies (see reviews 

; bvlSellwood fc Wi lkins(mlll993l: iKormendv fc Kennicutti 

! 120041: iGad otti 2009a|, and references therein). One of 
the expectations that emerges from both observation and 
theory, in the framework of secular evolution processes 
induced by bars, is the rejuvenation of the stellar popula- 
tion in the central structural component of disk galaxies. 
Bars are able to collect gas in the disk from within the 

[ bar ends to the central parts of the disk, supposedly help- 
ing the buildi ng of bulges through central star formation 

' episodes (see Athanassoula '2005") . Simulations such as 
those of [Athanassoula & Misiriotis (2002), but including 
gas, show that the transfer of gas to the center should be 
fast, « 10® yr (E. Athanassoula, priv. comm.). 

To date, there is evidence for an enhanced star for- 
mation rate, i.e. current star-forming activity, in the 
centers of barred galaxies, mostly from studi es of nuclear 
Hii regions. For instance, iHo et all (|1997| ) found that 
Ha emission line luminosities and equivalent widths are 
enhanced in barred galaxies, as compared to unbarred 
galaxies, when one c onsiders early- type di sk galaxies only 
(see also iHuang et al. 1996; Alonso-Hcr rero fc KnapenI 
[200lt iJogee et al.ll2005l: lEUison et al.ii20lK 

Direct evidence supporting bulge building by bars 
from the ages of stars in bulges has proven to be 
much more elusive. Studies based on integrated col- 
ors have to deal with uncertainties that arise from the 



age-met allicity degeneracy and effe cts of dust extinction 
(see e.g- fGadotti fc dos Aniosll200lD . On the other hand, 
studies based on stellar spectral an alysis have been, to 
date, handicapped by poor statis tics (jPeletier et al.ll2007l : 
iPerez fc Sanchez-B lazquezjl2011f ). 

In this Letter, we use spectral analysis techniques on 
a large and well defined sample of barred and unbarred 
disk galaxies, drawn from SDSS, allowing us to compare 
mean stellar ages of bulges in barred and unbarred galax- 
ies with unprecedented statistical significance. We de- 
scribe the most relevant features of our data in the next 
section. Section [3] describes relevant details of our spec- 
tral analysis techniques, and results are shown in Sect, 
m We present our main conclusions in Sect. [5l 

2. SAMPLE 

The sample u sed he re is based on the one stud- 
ied in iGadottil (l2009b[F1 . which contains all galaxies 
in SDSS Data Re lease 2 with stellar m asses larger 
than 101°Mq (from lKauffmann et al.ll2003a[ ). at redshift 
0.02 < z < 0.07, and with axial ratio b/a > 0.9. These cri- 
teria provide a sample which is both representative and 
suitable for 2D bulge/disc/bar decomposition, as select- 
ing face-on galaxies minimizes dust and projection effects 
and eases the identification of bars. The reader is referred 
to that paper for a detailed discu ssion on selection ef- 
fects. Through 2D decomposition, IGadottil ()2009b[ ) pro- 
vides reliable structural parameters, such as the total 

^ See |http : //www ■ sc. eso ■ org/- dgadott i/buddaonsdss ■ html | 
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stellar mass, bulge stellar mass, bulge effective radius 
Te, bulge/total and bar structural parameters, which are 
used in this Letter. In short, bulge stellar masses were 
obtained from bulge luminosity and mass-to- light ratio in 
the i-band, the latter derived from the bulge g — i color. 
Galaxy total stellar masses were obtained by adding the 
masses of its components. 

To verify whether a galaxy is barred, typical bar 
signatures were searched for through inspection of the 
galaxy image, isophotal contours and a pixel-by-pixel ra- 
dial intensity pro file (as in Scannapicco ct al. 2010 and 
ISheth et aLl(2010l ). It should be noted that due to the 
limited spatial resolution of SDSS images, we miss most 
bars with semimajor axis shorter than '^2 — 3 kpc, which 
are found mainly in very late-t ype spirals (later than Sc; 
lElmegreen fc ElmegreenI [19851 ). 

We selected from the previous sample all disk galax- 
ies with bulges and filtered them in terms of bulge-to- 
total and signal-to-noise ratios (see fg]) so that our effec- 
tive sample has 251 barred and 324 unbarred galaxies, 
of which 187 are AGNs ac cording to the classification 
in IKauffmann et all (|2003bl ) . No galaxy shows emission 
lines with equivalent width larger than 20A and there- 
fore, according to the criterium in IPetersoiil ()1997l table 
1.1), aU AGNs are type 2. 

The galaxy stellar mass distribution is similar for 
barred and unbarred galaxies, with a difference with 
statistical significance of less than Icr, according to a 
Kolmogorov-Smirnov test (hereafter KS). This means 
that possible biases, which would be caused by differ- 
ent mass distributions, are absent in our results. For 
instance, this indicates that global star formation histo- 
ries are similar in the barred as in the unb arred galaxies 
in our sample (see 'Kau ffmann et al.l l2003a). 

The spectra, obtained from the SDSS database , have 
a spectral resolution of A/AA - 1800 |York et al.ll2000ll 
and an average S/N in the spectral region covered by 
the SDSS 5-band of ~ 21. The g-band spectral region is 
most relevant for our analysis, as it contains the spec- 
tral features most sensitive to the stellar parameters we 
aim at deriving, i.e. ages and metallicities. The spectra 
were brought to restframe and corr e cted f or Galactic ex- 
tinction as in iCid Fernandes et all (|2005l ). It should be 
noted that SDSS spectra are taken through a fi xed fiber 
size on sky of 3" , centered at the galaxy center. iGadottil 
(|2009b) shows that for most galaxies in his sample the 
light within the fiber diameter is emitted mainly by bulge 
stars. In this stu dy, we haye used the results from the 
decompositions in IGadottil (|2009bD to measure the disk 
contribution. We verify, through KS tests, if the distri- 
butions of such disk contamination inside the fiber are 
similar for barred and unbarred galaxies, avoiding pos- 
sible related biases. This is discussed on a case-by-case 
basis. 

3. ANALYSIS 

We use the spectrum fitting code StarlighiEI 
(jCid Fernandes et al.ll2005[ ) to compare, on a pixel-by- 
pixel basis, the SDSS spectra to stellar population (SP) 
models. 

In short, STARLIGHT mixes different computational 
techniques to fit an observed spectrum 0\ with a com- 

^ See |http : //www. starlight .ufsc .br| 
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Figure 1. Example of a Starlight fit. The galaxy identification 
corresponds to SDSS plate-mjd-fiber numbers. Observed spectrum 
and model fit are shown as black and red lines, respectively. Pix- 
els not considered in the spectral fitting (either emission lines or 
clipped pixels), are marked as gray points. The average absolute 
deviation A = (J]|(Ma - Ox)/Ox\)/Npi^eis of the fit is given. 
The bottom panel shows the residuals 0\ — M\ . 



bination of N^, simple stellar population (SSP) models. 
Extinction is modelled as due to foreground dust and 
line-of-sight stellar motions are modelled by a Gaussian 
distribution G centered at velocity and with disper- 
sion tj* . Both kinematical and SP parameters are derived 
during the fit and the best model spectrum is given by: 




A^A ® G'(u*,f7*) 



(1) 



where hj^x is the spectrum of the jth SSP normalized at 
Ao, rx is the reddening term, x is the population vec- 
tor whose components Xj (j=l,...,7V*) represent the frac- 
tional contribution of each SSP to the total synthetic 
fiux at Ao , Mxq is the synthetic flux at the normalization 
wavelength, G'(w* , ct* ) is the line-of-sight stellar velocity 
distribution, and €5 denotes the convolution operator. 

Known regions of emission lines in AGNs were masked 
for the whole sample, whether the galaxy is an AGN or 
not, to ensure a homogeneous analysis. As only type 2 
AGN are present in our sample, the non-stellar AGN con- 
tribution is limited to a few percent of the flux and do not 
affect the SP parameters derived (|Cid Fernandes et al.l 
12004 Cid-Fernandes, priv. comm.). Bad pixels or 
sky background residuals are clipped during the fit, 
when pixels deviate by more than three times the 
r.m.s. between Ox and M x- We refer the reader to 
iCid Fernandes et al.l ()2005t ) and references therein for 
more details on Starlight. 
The SP models adopted are those of 

a 



range 3540-7400A at a resolution of FWHM ~ 2.5A. 
Ages range between 63Myr and 18 Gyr, and metallicities 
[Af/H] between -2.32 and -f 0.22. A random fit is shown 
in Fig.[l] where the observed spectra is shown in black 
and the Starlight model in red. 

4. RESULTS 

In Fig. [5] we show the results of our stellar popula- 
tion analysis, where light-weighted mean ages and metal- 
licites - normalized at wavelength 4020A - are plotted 
for bulges in barred and unbarred galaxies. 

The robus tness of Starlight results were thoroughly 
discussed in ICid Fernandes et al.l ()2005D . with tests per- 
formed also on SDSS spectra. They conclude that mea- 



iVazdekis et al.l (120101) with an updated stellar library 
((Falcon-B arroso et al.ll2011f) . and cover the wavelength 
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Figure 2. Ages vs. metallicities of bulges in unbarred and barred 
galaxies (top and bottom panels, respectively). Results for the 
whole sample are shown in the left-hand column, and for different 
bulge sizes in the remaining columns, as indicated. 



surement uncertainties - which can be parametrized by 
S/N - dominate the final errors. From their Table 1, 
the uncertainty on the light-weighted parameters are of 
~20% for a S/N= 10 spectrum. Although we have pro- 
cessed all spectra, we rely our analysis solely on the re- 
sults from spectra with S/N> 10, which constitute 93% 
of the sample. 

To circumvent our inability of detecting the short bars 
in galaxies with Hubble types later Sc, we remove from 
the analysis all galaxies with bulge-to-total luminosities 
ratios below 0.043, which is the typica l ratio for these 
latest Hubble types (jGraham fc Worlevl[2008.) . 

In Fig. [5] we show the results for the whole sample in 
the left-hand column, and separated into bulges with 
re < 0.8", 0.8" < re < 1.3", and r^ > 1.3" ii i the re- 
ma ining columns (re in the «-band, as given in iGadottil 
I20d9b l. We note that the parameter space covered by 
the results does not change significantly with bulges size, 
which is evidence that there are no significant biases from 
disk contamination inside the SDSS fiber, even for small 
bulges. Note that bulge light dominates over disk light 
through ss 2 ti mes the bulge r^ from the galaxy center , 
on average Csee lMoriondo et al.lE998l : lMorelh et al.ll2008l : 
iGadotti 2009bl). There seems to be a lack of very old 
(> 10 Gyr) bulges in unbarred galaxies with rg < 0.8". 
This is not the case for barred galaxies. If this was a con- 
sequence from disk contamination, one would see bulges 
in unbarred galaxies to have younger ages on average 
than in barred galaxies (opposite to what we find - see 
below). Moreover, stellar populations younger than ~ 5 
Gyr are seen equally in small and large bulges, barred 
and unbarred galaxies. It is for these ages that we find 
a significant difference between bulges in barred and un- 
barred galaxies. We thus conclude that our results are 
not affected by disk contamination. 

In Fig.|3] we present the normalized distributions of 
ages and metallicities derived, for an upper limit of bulge 
stellar mass of lO^'^'^^M© (above which there are no 
barred galaxies in our sample). The top panels show the 
distributions for the whole sample, and the middle and 
bottom panels show the distributions for AGNs and nor- 
mal galaxies, respectively. The age distribution of bulges 
in barred galaxies shows an excess of populations younger 
than '^4 Gyr. This feature is enhanced when we divide 
the sample in normal (non-active) galaxies and AGNs: 
the excess of young populations is better seen in the dis- 
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Figure 3. Normalized distributions of ages (left-hand column) 
and metallicities (righ-hand column) for bulges in barred and un- 
barred galaxies (red and black lines, respectively). Distributions 
for the whole sample, AGNs and non-active galaxies are given sep- 
arately. This sample has 251 barred and 324 unbarred galaxies 
(106 and 81 AGNs, respectively). 
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Figure 4. Normalized distributions of ages for bulges with 
logMjjujge < 10.1 on the left-hand side, and 10.1 < logA/jj^jge < 
10.85 on the right-hand side. The low-mass sub-sample has 160 
barred and 157 unbarred galaxies (56 and 25 AGNs, respectively) 
and the high-mass sub-sample has 91 barred and 167 unbarred 
galaxies (50 and 56 AGNs, respectively). 



tribution of normal galaxies, and disappears in AGNs. 
In fact, a KS test shows that the probabilities that the 
bulge mean age distributions are drawn from different 
populations for barred and unbarred galaxies are 98.65% 
for the whole sample, 50.13% for AGNs and 99.94% for 
normal galaxies only. Thus, the difference between the 



4 



P. Coelho & D. A. Gadotti 



mean stellar ages of bulges in non-active barred and un- 
barred galaxies is significant to a level of almost 4(t. 

The distributions of bulge stellar mass for the sample 
in Fig. |3] are, however, statistically different for barred 
and unbarred galaxies. Barred galaxies have less massive 
bulges than unbarred galaxies, even though their total 
stellar mass is similar, and also have larger disk con- 
tamination within the fiber. Therefore, and because less 
massive bulges tend to have younger mean ages, one can- 
not tell from Fig. [3]alone that bars do indeed turn bulges 
younger. 

For this reason, we inspected several hundred bulge 
mass intervals in search for those where the distributions 
of bulge mass, and the distributions of disk-to-total light 
ratios inside the fiber, are the same for barred and un- 
barred galaxies at ~ la level (ensured with KS tests). 
We could not find intervals where both AGNs and nor- 
mal galaxies show equal bulge mass distributions, and 
then focused on choosing optimal distributions for nor- 
mal galaxies only (a detailed study of the AGNs will 
follow in a separate paper). We thus came to two inter- 
vals of bulge stellar mass - below 10"'^'''^Mq and between 
lO^^'^M© and lO^^'^^M© - whose corresponding results 
are shown in Fig. |4l 

In the lower mass bin, the distributions of bulge ages 
are statistically similar, in contrast with Fig. [31 which 
refers to the whole sample. However, for the high-mass 
bin, the distribution of bulge ages for barred galaxies is 
clearly bimoda l. The mixture modelling statistical KMM 
test (see Ashman et al.l[l994f ) indicates that this is so at 
a confidence level of 99.9993%, i.e. > 4ct. The same test 
results in two normal distributions, with peaks at 4.7 
and 10.4 Gyr. We have run the KMM test in all other 
distributions discussed, and none resulted in a statistical 
significance larger than ~ Itr, in particular the distribu- 
tion of bulge age for unbarred galaxies. This bimodality 
does not seem to be a result from biased samples, and 
is a clear evidence of difference between the mean stellar 
ages in bulges of barred galaxies, as compared to un- 
barred galaxies. 

To better inspect this bimodality, we show in Fig. [5] 
the age distributions for several bulge mass intervals. 
This is essentially the same as having matched distri- 
butions of bulge mass, since the mass intervals at each 
panel are narrow. A first signal of bimodality in the age 
distribution of non-active galaxies appears in the interval 
log Mbuige between 9.7 and 10.2 (the KMM test yields a 
significance of 99.8928%) and the peaks reach compara- 
ble strengths, as in Fig. U (significance of 99.9813%), in 
the interval log Mbuige = 9.9 - 10.4. A similar analysis 
was done for AGNs, corroborating the results above that 
those effects are inexistent in AGNs. 

This result provides statistically-based corroboration 
that bars have important effects on the processes of bulge 
building, though the dependence of these effects with 
bulge mass, and the bimodality in the bulge age distri- 
bution, have yet to be explained by theoretical work. 

Furthermore, bars also help building a reservoir of fuel 
for AGN activity. In the low bulge mass bin, 35% of 
barred galaxies are AGN, whereas this fraction drops to 
only 16% when one considers unbarred galaxies. In the 
high mass bin, 55% of barred galaxies are AGN, whereas 
34-% of unbarred galaxies are AGN. In both mass bins, if 
a galaxy is barred it has a higher chance of hosting an 
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Figure 5. Normalized distributions for bulge ages in non-active 
galaxies, for several mass intervals as indicated. Distributions for 
barred and unbarred galaxies are sliown in red and black lines, 
respectively. 



AGN. In the high mass bin, processes such as mergers 
might be acting to help fueling AGN, diminishing the 
difference in the fraction of AGNs between barred and 
unbarred galaxies. Bars are not a necessary nor suffi- 
cient condition for a galaxy to host an AGN, but these 
results indicate that in some cases bars do help in fuel- 
ing AGN activity. There are several studies in the lit- 
erature - with opposing results - on this issue (see e.g. 
Mulchaev k Re gan 1991 iHo et al.l [I997I: iKnapen et all 
2000tlLaine et al.ji2002i) . but they agree that homogene- 
ity on the detection of AGN activity and the assessment 
on the presence of a bar, as well as sample selection, are 
critical. Our sample is carefully drawn from a volume- 
limited sample, the data set is homogeneous, and both 
the AGN and bar classifications are done throughout in 
a consistent fashion. Nevertheless, we underline that our 
sample comprises only massive galaxies and that short 
bars are mostly missed, as discussed above. 

Further work is necessary to better understand why 
the difference between the ages of bulges in barred and 
unbarred galaxies disappears in AGNs, and why the di- 
chotomy in the bulge ages of massive barred galaxies is 
not evident in AGNs. A likely interpretation is that feed- 
back from the AGN activity will at some point push gas 
back from the galaxy cen ter, preventing iiew st ar for- 
mation episodes (see e.g. iSchawinski et al.ll2007[) . We 
stress that even if bulges in barred AGNs are typically 
less massive, they are not younger than their unbarred 
counterparts. 

The metallicity distributions show no important dif- 
ference between barred and unbarred galaxies, but we 
will further explore the chemical enrichment in terms of 
a-elements over iron abundances ratios as a function of 
bulge morphology in a separate paper (P. Coelho & D. 
A. Gadotti, in preparation). 
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We have compared the ages with bar structural pa- 
rameters - effective surface brightness, effective radius, 
ellipticity, Sersic index, semi-major axis, boxiness and 
bar-to-total ratio - but found no evidence for a relation 
between any of the bar properties and the age of the bulge 
population. This is not necessarily surprising, as bulge 
building by bars depends on complex physical processes 
and time scales, and the availability of gas. For instance, 
a bar which was once strong, but is now weakened for 
any reason [bar weakening is more likely to happen than 
bar destruction, (e.g. [Athanassoula ct al. 2005)], could 
have contributed substantially to build a young popula- 
tion in the bulge of its host galaxy (if enough gas was 
available), and would now weaken a correlation between 
bar strength (e.g. ellipticity) and bulge age. 



5. CONCLUSIONS 

We derived stellar ages for a sample of 575 bulges in 
disk galaxies, 251 of those containing bars. When we con- 
sider the sample with bulges stellar masses < 1O^°'*^M0, 
we find that the mean stellar ages of bulges of barred 
galaxies are on average lower than that of unbarred 
galaxies, at a statistical significance of 99.94% (or al- 
most 4cr), when one considers non-active galaxies only. 
In this sample the galaxy mass distributions are similar 
between barred and unbarred galaxies. 

To make the distributions of bulge stellar mass of 
barred and unbarred galaxies similar, we split the sample 
in two bulge mass bins. We find that bulges in massive 
non-AGN barred galaxies (logM > lO^^'^M©) show a 
bimodal stellar age distribution, at a confidence lever of 
99.9993%, or more than 4cr. This can be described as 
two normal distributions, centered at 4.7 and 10.4 Gyr. 
This bimodality is present above a characteristic mass 
logMbuige = 9.7 — 10.2 and is absent for unbarred galax- 
ies or AGNs. As discussed above, this is a strong ob- 
servational evidence that corroborates scenarios of bulge 
building by secular evolution processes induced by bars. 
On the other hand, this bimodality, and the ages for the 
two distributions it consists of, are new constraints yet 
to be explained by successful theories of bar evolution. 

We have verified that our results are not caused by 
biases in the samples compared. We have taken into ac- 
count, at separate instances, the galaxy mass and bulge 
mass distributions of barred and unbarred galaxies, in or- 
der to have samples with e.g. similar star formation his- 
tories. We have also used sub-samples in which the con- 
tributions from disk light within the SDSS fiber, through 
which the spectral information used here is taken, are 
similarly distributed. Finally, we have not considered 
very late-type galaxies, for which the presence of a bar 
cannot be reliably assessed in our sample. Therefore, our 
results cannot be attributed to a bias in sample selection 
or a fiaw in the methodology we use. 

Finally, let us point out again that samples of barred 
and unbarred galaxies with similar galaxy total stellar 
mass distributions have statistically significantly differ- 
ent distributions of bulge mass, in the sense that bulges in 
barred galaxies have lower masses. This seems to be un- 



expected because one expects bulge building by bars and, 
secondly, bars grow from disks and thus one would ex- 
pect changes in the sense of lower disk masses. Progress 
in both theoretical and observational results is needed to 
clarify this trend. 
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